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Abstract 

Experimental observations of animal collective behavior have shown 
stunning evidence for the emergence of large-scale cooperative phe- 
nomena resembling phase transitions in physical systems. Indeed, 
quantitative studies have found scale-free correlations and critical be- 
havior consistent with the occurrence of continuous, second-order phase 
transitions. The Standard Vicsek Model (SVM), a minimal model of 
self-propelled particles in which their tendency to align with each other 
competes with perturbations controlled by a noise term, appears to 
capture the essential ingredients of critical flocking phenomena. In 
this paper, we review recent finite-size scaling and dynamical studies 
of the SVM, which present a full characterization of the continuous 
phase transition through dynamical and critical exponents. We also 
present a complex network analysis of SVM flocks and discuss the 
onset of ordering in connection with XY-like spin models. 



1 Introduction 

Nature offers abundant manifestations of collective motion phenomena in 
self-propelled living systems at all scales, from biomolecular micromotors, 
migrating cells and growing bacteria colonies, to insect swarms, fish schools, 
bird flocks, mammal herds and even human crowds. Instead of focusing on 
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the specific details that make each of these biological systems unique, statisti- 
cal physicists have been studying the general patterns of biological collective 
motion, aiming to identify the general laws and underlying principles that 
may govern their behavior pQ. 

From this perspective, one important question to address is the onset 
of ordered macroscopic phases, i.e. the way in which individuals having 
short-range interactions are capable of self-organizing into large-scale coop- 
erative patterns in the absence of leaders or other ordering cues from the 
environment. By analogy with large molecular systems, flocking and swarm- 
ing phenomena can be associated with phase transitions that depend on a 
few parameters that characterize the macroscopic states, such as the density 
of individuals and the flock size. For instance, Buhl et al. [2] investigated 
the collective motion of locusts, which display a density- driven transition 
from disordered movement of individuals within the group to highly aligned 
collective motion. In analogy with second-order (or continuous) phase tran- 
sitions in physical phenomena such as critical opalescence, in which liquids 
appear opaque due to density fluctuations at all wavelengths of visible light, 
the onset of a coordinated alignment in the motion of marching locusts was 
found to display second-order behavior. Moreover, a critical locust density 
for the occurrence of the continuous phase transition was identified. 

Familiar to most of us is also the rich variety of collective motion phenom- 
ena displayed by flocks of birds, where highly correlated long-range ordering 
effects are readily apparent (see e.g. [3] for some stunning footage of complex 
flocking patterns and other kinds of swarm behavior). Indeed, some stud- 
ies have quantitatively characterized the formation of ordered bird flocks 
in terms of continuous phase transitions. For instance, Cavagna et al. [I] 
obtained high resolution spatial data of thousands of starlings using stereo 
imaging in order to calculate the response of a large flock to external per- 
turbations. This study showed that behavioral correlations are scale-free, 
a signature of critical states in second-order phase transitions. Correlation 
functions revealing the pairwise interaction between homing pigeons have re- 
cently been investigated by means of high-resolution lightweight GPS devices 
as well [5J. 

On the theoretical side, Vicsek et al. |6] proposed a minimal model to 
study the onset of order in systems of self-driven individuals, which was later 
followed by other investigations by means of agent-based modeling [7J |8] , the 
Newtonian force-equation approach [HI [10], and the hydrodynamic approxi- 
mation [TTjEO]- This so-called Standard Vicsek Model (SVM), assumes that 
neighboring individuals tend to align their direction of movement when they 
are placed within a certain interaction range. This alignment rule, which 
would trivially lead to fully ordered collective motion, is complemented by a 
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second one that introduces noise in the communications among individuals. 
The SVM assumes that each individual can assess exactly the direction of its 
neighbors, but it incurs an error when adjusting its own direction of motion 
to match that of its neighbors' average. This type of noise, sometimes identi- 
fied as Angular Noise, was shown to drive the system through a second-order 
transition between the ordered phase of collective motion and a disordered 
phase. Although several variations to the SVM have later been considered in 
the literature, such as different noise types, models without alignment rule, 
adhesion between neighbor individuals, bipolar particles, etc (see PQ for a re- 
view), here we will chiefly focus on the original SVM as formulated by Vicsek 
et al. in their seminal paper [6]. 

In this context, the aim of this work is to review recent findings for the 
SVM that confirm the continuous nature of the order- disorder phase tran- 
sition driven by the noise amplitude, as well as to present a novel complex 
network analysis that provides new insight into the onset of ordering. The 
rest of this paper is organized as follows. In Section 2, we present the defini- 
tion of the SVM. Section 3 reviews recent finite-size scaling and dynamical 
studies of the SVM, leading to a fully quantitative characterization of the 
continuous phase transition through dynamical and critical exponents. In 
Section 4, we present a complex network analysis of SVM flocks, we study 
their topology, and we discuss the onset of ordering in connection with XY- 
like spin models. Finally, our Conclusions appear in Section 5. 

2 The Standard Vicsek Model (SVM) 

The SVM [S] consists of a fixed number of interacting particles, N, which are 
moving on a plane. In computer simulations, that plane is represented by 
a square of side L with periodic boundary conditions [HI [13] . The particles 
move off-lattice with constant and common speed vq = \v\. Each particle 
interacts locally adopting the direction of motion of the subsystem of neigh- 
boring particles (within an interaction circle of radius Rq centered in the 
considered particle), which is then perturbed by the presence of noise. Since 
the interaction radius is the same for all particles, we adopt the interaction 
radius as the unit of length throughout, i.e. R = 1. 

The evaluation of the average angle of motion of each individual's neigh- 
bors at time t, is affected by a noise term. Hence, the updated direction 
of motion for the i—th particle, is given by 

+ V&, (1) 
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where rj is the noise amplitude, the summation is carried over all particles 
within the interaction circle centered at the i— th particle, and is a realiza- 
tion of a 5-correlated white noise uniformly distributed between — tt and ir. 
The noise term can be thought of as due to the error committed by the indi- 
vidual when trying to adjust its direction of motion to the averaged direction 
of motion of its neighbors. 

The model dynamics is usually implemented by adopting the so-called 
backward update rule: after the position and orientation (velocity) of all 
particles are determined at time t, the position of the particles at time t+1 
is updated according to 

%i Xi -\- Vi , (2) 

which is then followed by the update of all velocities at time t + 1 according to 
Eq.flTJ. However, different updating schemes implemented in the literature 
may lead to spurious effects in the simulations [T4"] . 

The parameters of the model are the number of particles N, the dimension 
of the displacement space d, the linear size of the displacement space L, the 
interaction radius R , the particle density p, the particle speed v , and the 
noise amplitude rj. Notice, however, that p = N/L d , where d — 2 in the 
standard case, and that Ro = 1 can be chosen as the unit of length. Moreover, 
Vq 1 merely plays the role of a "thermalization parameter" that measures how 
many times, on average, two neighbors check out each other's positions while 
they remain at a distance within the unitary interaction radius. Therefore, 
relative large values of speed (typically v > 0.3) correspond to the low- 
thermalization regime, characterized by highly anisotropic diffusion and the 
manifestation of simulation artifacts in the form of directionally quantized 
density waves (see Ref. [15] for details). The relevant region of parameter 
space addressed in this paper is the high-thermalization regime for Vq < 0.1 
that guarantees isotropic diffusion. Within this regime, the choice of vq has 
only mild effects on the order parameter probability distributions and it does 
not affect the nature of the phase transition [T3] . 

3 Critical Behavior of the SVM: Finite-Size 
and Dynamic Scaling 

The SVM exhibits a far-from-equilibrium continuous phase transition be- 
tween ordered states of motion at low noise levels and disordered motion at 
high noise levels. The natural order parameter is given by the absolute value 
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of the normalized mean velocity of the system [HI ITS] : 



N 

1 , - 



1 N 

dE4 ( 3 ) 



where is close to zero in the disordered phase and close to unity in the 
ordered phase. In order to anlyze the critical nature of the phase transition 
exhibited by the SVM, we have studied behavior of the order parameter ip 
and other quantities related to the probability distribution of (p, such as the 
susceptibility (% = ^-[(y? 2 ) — (<^) 2 ]) and the Binder cumulant (U — 1— J^ 2 u )- 
The phase transition has been investigated by means of two independent 
approaches, namely a finite-size scaling analysis and a short-time critical 
dynamic technique [HI [T7] . 

Let us briefly mention that the second-order nature of the SVM far- 
from-equilibrium phase transition has been a matter of some debate. In 
fact, early simulations by Vicsek's group, which were consistent with a con- 
tinuous second-order transition [6J, were later challenged by Gregoire and 
Chate [IB] , who claimed that the transition was discontinuous, i.e. of first 
order. The controversy was further stimulated by subsequent papers of Vic- 
sek et al. [331 EE!, Aldana et al. 122 EJJ [22], Dossetti et al. [23], and 
Baglietto and Albano [131 HH EE], supporting the critical nature of the tran- 
sition, which were in conflict with additional results published by Chate and 
collaborators [2"lj [221 EE]- In the remainder of this work, we assume that 
the SVM (using the so-called Angular Noise type) undergoes a second-order 
phase transition. For further details, see Refs. [HI US [22] . 

While critical systems are characterized by the divergence of correlation 
length and time scales, computational simulations are limited to systems of 
finite extent and finite time observation intervals. In systems of finite size L 

(where, in our case, L = \J~^) the correlation length cannot grow beyond L. 

Instead of the theoretically expected power-laws at criticality (e.g. (p ~ \t\@ 
and x ~ l r l~ 7 ) where r represents the distance to the critical point and [3 and 
7 are critical exponents), for finite-size samples of critical systems, stationary 
quantities must obey scaling laws of the form: 

<p(T,L)=L~^(TL^), (4) 

X {t,L) = ISI v x{tIM v ), (5) 

where <p> and x are scaling functions, and v is the correlation length critical 
exponent. In the case of the SVM, a very good agreement between simulation 
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Figure 1: Plots of the order parameter ( (p) versus the noise amplitude (rj). 
Results obtained within the stationary regime for samples of different size 
and by varying the number of individuals, as listed in the figure [16J. 



results and Eqs.(j4])-([5]) are obtained. Fig. [T] shows the behavior of (p vs rj for 
different sizes and densities, while Fig. [2] displays the collapse obtained by 
applying Eq.pj to a known scaling law for the densities [16J. The values of 
j3 and v used in the simulations are listed in Table El 

On the other hand, when a critical system, initially placed either in a 
ground state or in the disordered phase, is suddenly brought to the critical 
point, it is also possible to measure its critical dynamic evolution. In this 
way, one can obtain not only the static exponents u, and 7) but also the 
dynamic exponent z that governs the time-dependent growth of the correla- 
tion length [17J. Table [3] summarizes the exponents obtained for the SVM 
by means of stationary and dynamic measurements [16J. 



Considering statistical errors, the obtained exponents satisfy the hyper- 
scaling relationship [dv — 2/3 ===== 7), which is firmly established for critical 
systems in equilibrium but less documented for the case of nonequilibrium 
critical systems (see e.g. |27]). This finding is a strong indication that the 
SVM undergoes a continuous nonequilibrium phase transition. It worth men- 
tioning that the set of critical exponents found for the SVM constitutes a new 
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Figure 2: Finite-size scaling analysis of data obtained under stationary condi- 
tions for different densities and numbers of individuals, showing log-log plots 
of the rescaled order parameter <p* = Lp[V)L^l v versus the density-rescaled 
noise rf = M^h^Mi^V 1 ' [16] . Here, the leftmost (rightmost) points corre- 
spond to data measured closest to (farthest from) the critical point. 

universality class. Also, the agreement between results obtained from sta- 
tionary and dynamic measurements provides strong support to the SVM's 
critical nature. 



4 The Onset of Orientation Ordering in SVM 
Flocks 

The SVM, which describes a far-from-equilibrium phenomenon, has often 
been compared with another archetypical case in the study of critical phe- 
nomena but under equilibrium conditions, namely the XY model [281 129] , In 
the XY model, nearest-neighbor interacting spins may adopt any possible ori- 
entation depending on the strength of the interactions and the temperature. 
By considering the SVM as a model of interacting spins that can undergo 
displacements in the direction of the spin, the basic difference between both 
models is precisely due to those displacements. In fact, other relevant sym- 
metries for the study of phase transitions, such as the dimensionality of the 
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Table 1: Exponents of the SVM as obtained by using different methods. 
The acronyms S, RD, and DDC refer to stationary measurements, relaxation 
dynamics from ordered states, and dynamic measurements starting from dis- 
ordered configurations, respectively [16]. 
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space, the nature of the order parameter, and the range of the interactions, 
are the same in both models. On the other hand, it is well known that, 
according to the Mermin- Wagner theorem [30[ I3T]. the XY model, as well as 
all other equilibrium systems in 2D with short-range interactions and 0(2) 
symmetry defined on translationally-invariant substrates, cannot display or- 
dered phases. Therefore, the onset of ordering in the SVM in 2D is quite 
intriguing. What is, then, the essential ingredient that generates the ordered 
phase of SVM flocks, despite Mermin- Wagner? Naturally, the possible an- 
swers lie in the conditions required by the theorem to apply, which are not 
met by the SVM, namely: 

• the SVM is non-equilibrium (e.g. particle interactions do not conserve 
momentum) 

• SVM spins move and originate effective long-range interactions (flocks 
merge and dismember, carrying information across the system) 

• The SVM substrate is not translationally invariant (particle positions 
are distributed inhomogeneously in space) 

It is worth mentioning that the XY model defined on ID small-world 
networks (generated by the addition of a small number of random long-range 
"shortcuts" on top of a regular ID lattice) was found capable of sustaining 
ordered phases [32j [33]. Hence, the inhomogeneous arrangement of parti- 
cle positions in the 2D displacement space, interpreted as a complex net- 
work substrate in which "spins" undergo ferromagnetic interactions, could 
arguably be the key ingredient that allows the onset of ordered clusters in 
the SVM. 
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In order to explore this further, let us define the concept of flocks as com- 
plex network structures. Starting with a disordered initial state in which indi- 
viduals are randomly located within the LxL displacement space, SVM rules 
lead to the formation of local structures of interacting individuals. These 
structures, which we call flocks or clusters, are not permanent: their shape 
and size evolve with time, with new individuals and sub-flocks merging with 
them and, conversely, other individuals and sub-flocks separating from them. 
Moreover, besides the statistical perspective in which flocks are the funda- 
mental building blocks of the ensembles that characterize SVM stationary 
states, flocks can be regarded as "domains" that carry key information on 
the ordering of the system at the mesoscopic level. 

In order to gain insight into the structure of the clusters, let us evaluate 
the average path length, APL [31]. According to the standard definition 
used in the study of complex networks, for each pair of nodes or individuals 
(A, B) belonging to the same cluster, the path length Iab (also known as 
chemical distance) is given by the minimum number of links that one has 
to use in order to pass from one node to the other. By calculating all the 
pairwise node-to-node path lengths in the cluster and taking the average, one 
obtains the APL, which consequently is a characteristic length of the cluster. 
In Euclidean lattices, the volume of an object is related to its characteristic 
length by an integer power, i.e. the dimension of the object. Based on 
this observation, as well as on the experience gained in the study of fractal 
objects, it is customary to define the dimension (D) of a complex network 
according to: 

APL oc m x J D , (6) 

where m c is the complex network size or, in the present context, the cluster 
mass [31]. Figure [3J shows a log-log plot of the APL versus m c for SVM 
flocks corresponding to the critical point of the second-order phase transition. 
The best fit to the data yields D = 4.0(2), which strongly suggests that 
the effective dimension of the clusters is D — 4. Other complex network 
measures, not shown here for the sake of space, confirm that SVM flocks are 
characterized by a 4— dimensional topology [35] . 

As mentioned above, flocks can be regarded as domains that carry im- 
portant information on the ordering of the system at the mesoscopic level. 
Therefore, let us now study the onset of order within individual flocks as a 
necessary condition to have macroscopic ordering system-wide. Notice, how- 
ever, that the detailed mechanisms leading to the emergence of global order 
from locally ordered clusters are not well understood yet and remain an open 
question that lies beyond the scope of the present analysis. 
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Figure 3: Log-log plot of the average path length APL as a function of the 
cluster size m c at criticality. Fits of Eq.flH]) to the data yield D = 4.0(2). 

With this aim, here we investigate whether the topology of so-called frozen 
clusters, once particle displacements and cluster rearrangements are sup- 
pressed, is capable by itself of supporting the existence of an orient at ionally 
ordered phase. In our approach, we first generate configurations of clusters 
by applying the full SVM dynamics. Once the non-equilibrium stationary 
state is reached, we identify the clusters and "freeze" them, i.e. we disal- 
low any further displacements of the individuals. From that point on, the 
orientation of the particles is allowed to evolve according to Eq.flTJ, but sub- 
sequent displacements (Eq.©) do not occur. We will refer to this stage as 
"restricted SVM dynamics" . Since the full SVM dynamics has an entangle- 
ment between particle displacements and XY-type interactions, by resorting 
to "frozen clusters" we disentangle these two major components. This proce- 
dure allows us to explore the relation between local topology and the ability 
for the restricted SVM dynamics to sustain local ordering. 

Figure H] shows the dependence of the order parameter f on the inverse 
cluster mass m~ l . The clusters were first generated using the full SVM dy- 
namics with critical noise. After freezing them, the restricted SVM dynamics 
was applied using different noise values in the < rjf < 1 range, as indicated. 
Even for very large noise amplitudes, the order parameter in the large clus- 
ter limit (m c — > oo) tends to finite values, e.g. ip ~ 0.04 for rif = 0.9. 
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Figure 4: Log-log plot of the order parameter as a function of the inverse 
cluster mass for frozen clusters with different noise levels, as indicated. Inset: 
Plot of the asymptotic values of the order parameter, yjoo, versus the noise 
amplitude rif. The solid line shows the mean field (i.e. fully-connected graph) 
results (fMF = s i n ( 7rr ?/)/ 7r7 ?/ [23]. 



Although we are not presenting here a full finite-size scaling analysis, these 
results were obtained for a large system size (N = 4 x 10 4 ) and they provide 
strong evidence that frozen clusters are capable of sustaining order (<p > 0) 
in the thermodynamic limit (N — > oo). The Inset to Figure H] shows the 
order parameter extrapolations to the m c — > oo limit, (p^, as a function of 
the noise amplitude rif. For the sake of comparison, the Inset also shows 
exact results from the mean field solution obtained for an infinite density of 
individuals [23J, which closely follow the trend of our computer simulation 
results. Recalling that the effective dimension of SVM clusters is D = 4, we 
argue that the observed mean-field-like behavior is related to the fact that 
D = 4 is the upper-critical dimension of the XY model, which essentially has 
the same symmetries as the SVM defined on frozen clusters. 

Therefore, our findings indicate that the full SVM dynamics creates flocks 
with enhanced connectivity, which behave as 4-dimensional objects compact- 
ified into a 2-dimensional displacement space (more details will be published 
elsewhere [22] )• By virtue of this higher-dimensional complex network sub- 
strate, the approximately equivalent XY representation (attained by our 
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"frozen cluster" procedure) displays mean-field-like behavior, thus solving 
the apparent paradox posed by the Mermin- Wagner theorem. 

5 Conclusions 

In this work, we revisit the Standard Vicsek Model (SVM), a minimal model 
of self-propelled particles that displays a noise-driven phase transition be- 
tween the ordered phase (in which individuals align with each other and move 
coherently system-wide) and the disordered phase (in which the net trans- 
lational motion averages out to zero). By means of a combined finite-size 
scaling and dynamical analysis of the SVM, we present a full characteriza- 
tion through dynamical and critical exponents and confirm the nature of the 
SVM phase transition as continuous (second-order), in agreement with the 
early results by Vicsek et al. [6]. Finally, we analyze the onset of ordering in 
SVM flocks by performing a complex network analysis and interpreting the 
results in connection with XY-like spin models. 

In a forthcoming paper [33], we will extend the complex network anal- 
ysis of Vicsek flocks to variants of the model, in which perturbations are 
implemented either as Angular Noise (AN) or as Vectorial Noise (VN). We 
believe that this topological approach will provide key insight into the dif- 
ferent nature of the ordering phase transitions observed in the Vicsek model 
under the two noise types, namely continuous second-order for Angular Noise 
and first-order for Vectorial Noise. Indeed, the connection between the order 
of the transition and the topology of the flocks remains an important open 
question that certainly deserves further investigation. 
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